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Without the use of any empirical fitting to experimental or high-level ab initio data, we present a 
double-hybrid density functional approximation for the exchange-correlation energy, combining the 
exact Hartree-Fock exchange and second-order M0ller-Plesset (MP2) correlation with the Perdew- 
Burke-Ernzerhof (PBE) functional. This functional, denoted as PBEO-2, is shown to be accurate for 
a wide range of applications, when compared with other functionals and the ab initio MP2 method. 
The qualitative failures of conventional density functional approximations, such as self-interaction 
error and noncovalent interaction error, are significantly reduced by PBEO-2. 



I. INTRODUCTION 

Over the past two decades, Kohn-Shani density func- 
tional theory (KS-DFT) [D] has become a very popular 
electronic structure method for molecular and solid-state 
systems, due to its computational efficiency and reason- 
able accuracy. Although the exact exchange-correlation 
(XC) density functional E^clp] in KS-DFT remains un- 
known, accurate approximations to Exc[p\ have been con- 
tinuously developed to extend the applicability of KS- 
DFT to a wide range of systems [1| . 

The hierarchy of approximations to Exc[p\ has been 
formulated as a Jacob's ladder connecting the earth 
(Hartree theory) to the heaven (chemical accuracy) 
The first rung of the ladder is the local density approx- 
imation (LDA) [1, , representing the XC energy den- 
sity by the local density. Going beyond the LDA, there 
has been an increasing interest in the "parameter-free" 
functionals developed by Perdew and co-workers, such as 
PBE [6] and TPSS j7], demonstrating the usefulness of 
these functionals. As LDA, PBE, and TPSS belong to 
the first, second, and third rungs of the ladder, respec- 
tively, they systematically improve upon the description 
of short-range XC effects by climbing up the ladder (with 
slightly increasing computational costs). However, due 
to the lack of accurate treatment of nonlocal XC effects, 
these semilocal functionals (first three rungs) can lead to 
significant errors, characterized by self-interaction error 
(SIE), noncovalent interaction error (NCIE), and static 
correlation error (SCE), in situations where these failures 
occur [3,i|. 

The SIEs of density functional approximations (DFAs) 
can be reduced by hybrid DFT methods (fourth rung) 
P, combining a fraction of the exact Hartree-Fock (HF) 
exchange with a semilocal functional. As the fraction 
of HF exchange included in a hybrid functional can be 
rationalized by perturbation theory arguments p^ . the 
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PBE-based hybrid functional PBEO [TT|, has gradually 
gained its popularity. In view of the successive evolution 
of the PBE-based functionals: LDA (first rung) PBE 
(second rung) -J> TPSS (third rung) PBEO (fourth 
rung) ? (fifth rung), naturally, the question is, what 
could be a functional on the fifth rung (highest rung) 
of the ladder? In this work, we propose a PBE-based 
double- hybrid functional (fifth rung), and demonstrate 
its superiority to the functionals on the first four rungs 
of the ladder, for a wide range of applications. 



II. THE PBEO-2 FUNCTIONAL 

The SIEs and NCIEs associated with DFAs can be si- 
multaneously reduced by double-hybrid methods (fifth 
rung) |12h17| , combining a fraction of HF exchange and 
a fraction of second-order M0ller-Plesset (MP2) corre- 
lation |18|] with a semilocal functional. Accordingly, a 
PBE-based double-hybrid functional is given by 
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where Ef'^ is the HF exchange, E^^^ is the PBE ex- 
change, E^^^ is the PBE correlation, and E^^"^ is the 
MP2 correlation (a perturbative term evaluated with the 
orbitals obtained using the first three terms of the en- 
ergy). In this work, the two mixing parameters, ax and 
flc, are determined by physical arguments (i.e., without 
the use of any empirical fitting to experimental or high- 
level ah initio data). 

Based on the double-hybrid approximations proposed 
by Sharkas, Toulouse, and Savin [l^, and the closely 
related work of Bremond and Adamo [1^, Toulouse et 
al. 21] have proposed the linearly scaled one-parameter 
double-hybrid (LSIDH) approximation (see Eq. (10) of 
Ref. plj), showing the relationship between the mixing 
parameters, ax and ac, in a double- hybrid functional, 
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Applying Eq. ([2]) to Eq. ([ij, only one of the mixing pa- 
rameters needs to be determined. 

Functionals on the first four rungs of the ladder fail to 
describe long-range van der Waals (vdW) interactions in 
regions where there is no overlap of the subsystems, due 
to the lack of long-range correlation effects jj, HH • This 
problem can be greatly reduced by the DFT-D (KS-DFT 
with empirical dispersion corrections) schemes pTi . [23l - 
[25I or by a fully nonlocal density functional for vdW in- 
teractions (vdW-DF) [26j, showing generally acceptable 
performance on several noncovalent systems [l^] ■ 

By contrast, the MP2 correlation in a double-hybrid 
functional is itself fully nonlocal, and thus responsible 
for describing long-range vdW interactions. To deter- 
mine the fraction Oc of MP2 correlation in Eq. ([T]) , con- 
sider the two limiting cases where Oc = {ax)^ = and 
Qc = {o,x)^ — 1, which reduce to PBE and MP2, respec- 
tively. In view of the underbinding tendency of PBE and 
overbinding tendency of MP 2 in describing the interac- 
tion energies of many vdW complexes, a mixing of 50% 
of PBE correlation with 50% of MP2 correlation seems 
physically justified (as vdW interactions are inherently 
correlation effects). Employing this half-and-half mixing 
of the PBE and MP2 correlation functionals with the 
LSIDH approximation in Eq. we obtain 

= (1/2)1/3 ^ 0.793701. (3) 

Our resulting PBE-based double-hybrid functional, de- 
fined in Eqs. H]) and ([3]), is denoted as PBEO-2, where 
the "0" refers to no fitting parameters, and the "-2" refers 
to the post-KS treatment for the MP2 correlation. 
By contrast, the mixing parameters = 1/2 and 
= (1/2)3 ^ i^g^ g^j.g adopted in PBEO-DH an- 
other PBE-based double-hybrid functional (also based 
on Eqs. ^ and wherein a half-and-half mixing 

of the PBE and HF exchange functionals is effectively 
employed. Therefore, a comparison between the perfor- 
mance of PBEO-2 and PBEO-DH is particularly impor- 
tant, as it shows the significance of the physical argu- 
ments used in determining the ax and Oc- 

III. RESULTS 

For a comprehensive comparison of different func- 
tionals, we examine the performance of PBEO-2, other 
PBE-based functionals (LDA [JL15i], PBE TPSS Q, 
PBEO [ll|, and PBEO-DH [2^), three popular func- 
tionals (B3LYP (a hybrid functional) [P,^, B2PLYP 
(a double-hybrid functional) [ll|, and B2PLYP-D3 (a 
double-hybrid functional with empirical dispersion cor- 
rections) [TtI), and the ah initio MP2 method [3, on 
various test sets involvin g th e 223 atomization energies 
(AEs) of the G3/99 set lli-llH, the 40 ionization po- 
tentials (IPs), 25 electron affinities (EAs), and 8 proton 
affinities (PAs) of the G2-1 set [sl], the 76 barrier heights 
of the NHTBH38/04 and HTBH38/04 sets 'SS',!!, the 
22 noncovalent interactions of the S22 set |35|, |36j], the 66 



noncovalent interactions of the S66 set two interac- 
tion energy curves of weakly bound complexes, isodesmic 
reaction energies, and two dissociation curves of symmet- 
ric radical cations, with a development version of Q-Chem 
3.2 [38]. Detailed information about some of the test sets 
can be found in Ref . 39] . The error for each entry is de- 
fined as (error = theoretical value — reference value). 
Spin-restricted theory is used for singlet states and spin- 
unrestricted theory for others. Results are computed us- 
ing the 6-311-1— |-G(3df,3pd) basis set, unless noted oth- 
erwise. For efficiency, the resolution-of-identity (RI) ap- 
proximation [13] is used for calculations with the MP2 
correlation (using sufficiently large auxiliary basis sets). 
The interaction energies of the weakly bound systems 
are computed with the counterpoise correction '411] for 
the basis set superposition errors (BSSE). 

As shown in Table [H for the AEs of the G3/99 set, 
B2PLYP and B2PLYP-D3 have the best performance 
(partly due to the fact that this test set overlaps with 
their training sets). PBEO-2 performs similarly to PBEO- 
DH, PBEO, TPSS, and B3LYP, and significantly out- 
performs LDA, PBE, and MP2. For the IPs, EAs, and 
PAs of the G2-1 set, LDA performs poorly, while all the 
other functionals perform similarly to MP2. For the bar- 
rier heights of the NHTBH38/04 and HTBH38/04 sets, 
PBEO-2, PBEO-DH, B2PLYP, and B2PLYP-D3 have 
comparable performance, and are the best functionals 
here [13] . For the noncovalent interactions of the S22 and 
S66 sets, B2PLYP-D3 performs the best (mainly due to 
the empirical dispersion corrections). PBEO-2 performs 
similarly to MP2, and significantly outperforms all the 
other functionals [4^. 

Here, we calculate the interaction energy curves for the 
parallel-displaced (in Fig. [T|) and T-shaped (in Fig. [5]) 
configurations of the benzene dimer as functions of the 
intermonomer distance R (defined in Ref. [37|V where 
the optimized geometries and reference values are taken 
from the S66x8 set [37]. Overall, the predicted interac- 
tion energy curves of B2PLYP-D3 and PBEO-2 are the 
closest ones to the reference curves (within an error of 1 
kcal/mol). LDA predicts too short equilibrium distances, 
and all the other functionals predict severely underbind- 
ing or even repulsive curves. MP2 incorrectly predicts the 
parallel-displaced configuration to be more stable than 
the T-shaped one. 

The isodesmic reaction energies of n-alkanes to ethane 
have been known to giv e sy stematic errors in conven- 
tional DFT calculations [43l447{ . The considered bond 
separation reaction, 

CH3(CH2)„CH3 + mCH4 ^ {m + 1)C2H6 (4) 

was first proposed by Redfern et al. fis*] to show the 
errors of B3LYP systematically increase as the number 
of chain units in n-alkanes, to, increases. To assess how 
PBEO-2 improves upon this problem, we calculate the 
isodesmic reaction energies of n-alkanes to ethane, where 
the optimized geometries and reference values are taken 
from Ref. [4f>]. As shown in Fig. [31 the predicted reac- 
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FIG. 1. Interaction energy curve for the parallel-displaced 
configuration of the benzene dimer as a function of the inter- 
monomer distance R (defined in Ref. [131), where Ro is the 
equilibrium distance. 



tion energies of PBEO-2 are extremely close to the refer- 
ence values (within an error of 0.02 kcal/mol), due to its 
accurate treatment of medium-range correlation effects 
(46j . By contrast, all the other functionals and MP2 
exhibit systematic errors with the increase of m. Note 
that B2PLYP-D3 overcorrects the systematic errors of 
B2PLYP. 

Due to the severe SIEs of semilocal functionals, spu- 
rious fractional charge dissociation can occur, especially 
for symmetric charged radicals [4^ , such as He J and 
ArJ. To examine how PBEO-2 improves upon the SIE 
problem, spin-unrestricted calculations, using the aug- 
cc-pVQZ basis set, are performed. The DFT results are 
compared with results from MP2 theory and the highly 
accurate CCSD(T) theory 49]. As can be seen in Fig. 
m the predicted HeJ binding energy curve of PBEO-2 is 
very close to the CCSD(T) curve. It appears that the SIE 
associated with PBEO-2 is more than three times smaller 
than the next best DFT method shown. For ArJ, PBEO- 
2 can dissociate it correctly (see Fig. [5]), which is a very 
promising result. However, a discontinuity undesirably 
appears in the derivative of PBEO-2 binding energy curve 
for ArJ, due to the use of MP2 correlation (as discussed 
in Refs. [SlI^I). 



FIG. 2. Same as Fig. [T] but for the T-shaped configuration. 



IV. CONCLUSIONS 

In summary, we have developed a PBE-based double- 
hybrid functional, employing a half-and-half mixing of 
the PBE and MP2 correlation functionals with the 
LSIDH approximation in Eq. Owing to its signif- 
icant improvement over LDA, PBE, TPSS, and PBEO, 
this functional, denominated PBEO-2, fits well into the 
fifth rung of the ladder. Our results indicate that PBEO-2 
is generally comparable or superior to PBEO-DH in per- 




FIG. 3. Errors for isodesmic reaction energies of n-alkanes to 
ethane. 
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FIG. 4. Dissociation curve of HeJ. Zero level is set to £^(He) 
+ _B(He^) for each method. 



4 



TABLE I. Mean absolute errors (in kcal/mol) of 10 methods for various test sets (see text for details). 



System 






PBE-based functionals 






Other functionals 


Ab initio 


LDA 


PBE 


TPSS 


PBED 


PBEO-DH 


FBEO-2 


B3LYP 


B2PLYP 


B2PLYP-D3 


MP2 


Cjj/yy (226) 


121.26 


21.51 


5.30 


6.29 


5.11 


5.82 


5.48 


3.97 


2.98 


9.18 


IP (40) 


5.73 


3.45 


3.59 


3.48 


3.17 


2.38 


3.79 


2.35 


2.35 


3.69 


EA (25) 


5.03 


2.41 


2.61 


3.10 


3.55 


3.01 


2.38 


2.16 


2.17 


3.95 


PA (8) 


5.91 


1.60 


1.81 


1.14 


1.09 


0.96 


1.16 


1.18 


1.13 


0.96 


NHTBH (38) 


12.62 


8.62 


9.16 


3.63 


1.57 


2.44 


4.69 


2.19 


2.44 


5.48 


HTBH (38) 


17.95 


9.67 


8.68 


4.60 


2.01 


1.39 


4.56 


2.17 


2.45 


3.38 


S22 (22) 


2.03 


2.72 


3.56 


2.46 


1.78 


0.61 


3.90 


2.25 


0.22 


0.89 


S66 (66) 


1.85 


2.23 


3.09 


2.10 


1.58 


0.67 


3.33 


1.96 


0.21 


0.45 




2 4 6 8 10 



R (angstrom) 
FIG. 5. Same as Fig.H but for Ar+. 



formance. As PBEO-2 contains a very large fraction (« 
79 %) of HF exchange and a large fraction (50 %) of MP2 
correlation, the SIE and NCIE associated with PBEO-2 
are significantly reduced, reflecting its superior perfor- 
mance for self-interaction problems and noncovalent in- 
teractions, respectively. However, due to the perturba- 
tive treatment of nonlocal correlation, the SCE of PBEO- 
2 can be enormous in situations where strong static cor- 
relation effects are pronounced [2, 8]. It remains chal- 
lenging to develop a generally accurate density functional 
resolving all the qualitative failures of semilocal approx- 
imations at a reasonable computational cost. 
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